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Time-dependent Ginzburg-Landau model with noise for neutral s-wave superconductor is derived 
from real-time finite-temperature quantum field theory by integrating over fermions in the density 
matrix propagator. Quantum decoherence due to environment of pair breaking fluctuations is iden- 
tified. Dynamics is described by Langevin equation in the near critical regime and by a nonlinear 
Schrodinger equation at zero temperature. 
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BACKGROUND 

Ever since the advent of the Ginzburg-Landau theory of superconductivity (l) there was continuous effort to derive 
it from the microscopic theory. The first step was made by Gorkov j| in the static case. The static case is of much 
thermodynamical interest, see e.g. (3| for recent results. However, dynamical out of equilibrium phenomena require 
a time-dependent version of the Ginzburg-Landau model. This include for example formation of topological defects 
during a fast second order phase transition ||| . Models of defect formation normally assume the dynamics is described 
in the near-critical regime by a classical Langevin-type partial differential equation with external white noise. First 
derivation of a time-dependent Ginzburg-Landau model was given in || , similar approach is described in the recent 
review Neither of these works explicitly discusses noise or decoherence of the order parameter field. This paper 
aims at filling this gap for the simplest example of neutral s-wave superconductor. 

A topic in its own right is an effective quantum action at zero temperature. Natural expectation is a kind of 
nonlinear Schrodinger equation analogous to the Gross-Pitaevskii equation for Bose-Einstein condensates. Such 
equations were indeed constructed for the phase field or in pcrturbative regime. An effective nonperturhative 
Lagrangian for the pair field A is one of our results. 

We construct the density matrix propagator for electrons with the help of the finite temperature Kyeldysh contour. 
The Hubbard-Stratonovich transformation is performed to introduce the order parameter A and at the same time to 
make electronic Hamiltonian quadratic. Next, we integrate over fermions and obtain effective description in terms 
of the auxiliary A field. From the effective action in terms of A's living on the time-ordered and anti-time-ordered 
branch of the contour one can read the effective equation satisfied by the order parameter including dissipation and 
noise. We identify decoherence due to interaction of the pair field with pair-breaking fermionic fluctuations. 

Our paper is organized as follows. In Section I we define the model and perform Hubbard-Stratonovich transfor- 
mation to introduce the pair field A. Our aim is to integrate out fermions in the density matrix propagator and read 
from it an effective action for A. High temperature case is considered in Section II. In this regime the expectation 
value of the order parameter is zero. In the spirit of the Ginzburg-Landau model we make an expansion in powers 
of A to obtain, in the long wavelength limit, partial differential Langevin equation with white noise. At the same 
time the long wavelength modes of A are subject to quantum decoherence and thus they can be considered as clas- 
sical. This equation can be extended below the critical temperature T c to describe spinodal decomposition as long 
as temperature remains greater than the developing gap. Further below T c the gap in the quasiparticle spectrum 
can no longer be neglected. In Section III wc perform an expansion around uniform equilibrium order parameter. 
The fluctuations decohere, their evolution is effectively classical. They satisfy integro-differential Langevin equation 
with colored noise and nonlocal dissipation kernel. Properties of the system change qualitatively as we approach zero 
temperature, see Section IV. Some sources of noise and decoherence are suppressed when temperature becomes much 
lower than the gap, other are restricted only to frequencies greater than twice the gap. Thus in the low temperature 
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and long wavelength limit there is no noise or decoherence, we obtain an effective quantum Lagrangian, which leads 
to a kind of nonlinear Schrodingcr equation for the order parameter A. 



I. THE MODEL 

S-wave neutral superconductor is described by the Lagrangian 



(1) 



ipyi are spin up/down electronic Grassman fields. V is a coupling constant for the s-wave pairing between spin-up 
and spin-down electrons. By the standard Hubbard-Stratonovich transformation the Lagrangian is equivalent to 



Lty, A] = f d 3 x (^Mc + ^-^VVa + rt&a) - ^A*A + (2tt) 3 / 2 A*^ t + (2tt) 3 / 2 A^{ 



(2) 



Functional integration over the auxiliary A leads back to L e . Variation of L with respect to the pair field A* gives 

V 



A 



(2tt) 3 /< 



(3) 



Note that our A differs from the standard definition by the factor (2tt) 3 / 2 . Our goal is to integrate over 
ip\il> to obtain effective theory in terms of pair fields. We split this task into subcases. 



II. ABOVE T c 

The equilibrium value of A vanishes above T c , we perform expansion around equilibrium in powers of A. It is 
convenient to introduce Fourier transforms in space 



A(t,x) = J d 3 k-^^A(t,k) . 



(4) 



(2^) 3 / 2 

For clarity we denote the transform of A by the same letter. In Fourier space the Lagrangian reads 

L[a,A]=L [a] + L A [A]+L int [a,A] , (5) 



^oH = J 



d 3 k 



^ (ia1d t a a - e^a) a a a ) 

«=T,I 





' A*A" 


j d 3 k 






V 



Hnt[a,A]= J d 3 k [A^M + AF^o] 



h 

where = ^ — fi and the interaction vertex i 



is 



F n [a] = J d 3 pa l (t,p)a ] (t,k-p) . 



(6) 



The density matrix propagator for A can be obtained by integrating over fermions. To this end we put the fields 
on the Kyeldysh contour. In the complex time plane the contour runs from to tf (time ordered branch), then 
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back from t/ to tj (anti-time ordered branch) and from f, to ^ — where (3 = \/ksT (imaginary branch). The 
imaginary branch prepares the system at ti in the state of thermal and chemical equilibrium characterized by p and 
/i. When tf/i is moved to ±00, we are left with fields living on the time-ordered branch (a + , A + ) and those on the 
anti-time-ordered branch (a_, A_). In terms of these the density matrix propagator can be written as 



DA* + DA + DA*_DA^e l I dt {LAlA+] ' L&lA - ]) e tS ' F[A+ ' A - ] 

e iS F [A+A-] _ J Da\_Da + Da* Da-e* I dt i L o[a+]+^int[a + ,A+]-L [a-]-L int [a-,A-]) 

In the absence of A the free fermionic operators of the same spin have propagators 

<f[ o+(t, k) al(t', k') ]), = S(k - fc>-^+^+^K*-0 [ 6( t 0/(-%_„J _ e(t > t)f<fl* % _J) 

(f [ a _(t, fc) ai(f , fc') = - ^e-^s+^+^K*-*') - *)/(-%_,„*) _ 0( t _ t')/(%_ mff ) 

(f [ o_(*, fc) at (*', k') }h = *(* - k')e^ + ^ + ^-^f(-Pe^ n J , 
(f [ a + (i,fc) aUt',k') }), = -Sik-k^e-^n+M+^-^ftfe^) . 



(7) 



(8) 

/(x) = (1 + e 31 ) -1 comes from the Fermi-Dirac distribution and T means ordering along the Kyeldysh contour. (. . 
means an average over the initial state which is a state of thermal equilibrium in a reference frame moving with 
velocity v. We take such a generalized initial distribution to keep trace of Galilean symmetry. 

The influence action Sp in the influence functional expiSp, compare Eq.(^), can not be found exactly because it is 
quadratic in the fermionic fields. Instead we work out Sp by perturbative expansion in powers of the order parameter 
A, A*. This expansion is in the spirit of the Ginzburg-Landau model. As (Fg[a])# = the leading order term is 
quadratic 



4 2) [A+,A_] 



dtdt' I d 3 kd 3 k' ( [2A+(t,k) (T[F£[a+(t)] Fp[a+(t')} ])# A*(t',fc') 

-A + (t,k) (f[ Ft[a+(t)] F p [a^t')] ]), A*_(t' , %') 
-Al(t,k) (f[ F- k [a + (t)} Ft [ a _(0] ])v A_(t',fc') ; 
+ [+ - -] ) • 



This expression can be worked out with the help of the correlators (||) to give 



? (2) 



A+,A_] 



2 / dtud 3 k A («, fc) A£(co, fc) G ab (w, fc) 
0,6=+,- 



(9) 



(10) 



C++ (cj — 2fcw — mo 2 , fc — 2m?7) 

G (lu — 2kv — bid 2 , k — 2mv) 

(u> — 2A;t7 — m» 2 , k — 2mz7) 

G \-{oj — 2kv — mi) 2 , k — 2mv) 



d 3 p 
d 3 p 



f(-(3e f )f(-f3e % _ ? ) /(^ /(/%_„.) 



u + £ p + H-p ~ tr ) u + e v + H-p + ir i 

f(pe ff )f((3e n _ f ) /(-/^/(-/^ 



w + € p + e k-p - lr l w + € p + e k-p + »»? 
2™ | d 3 p *( W + 6p + ejE _ p -) /(-#*)/(-/%-*) 
2m [ d 3 p S(cj + e p +er.) /(/^)/(/%_ J 



(11) 
(12) 



Where 77 is an infinitesimal positive parameter. G a b(uj, k) depends on v through ui + 2kv + and fc -I- 2mv which 
means that A transforms under Galilean symmetry like mass 2m pair field. We introduce an average and relative 
field, which are the usual coordinates used in Wigner transform, 



A, 



R 

A± 2 



(13) 



With this change of variables and with the identity (y + irf) 1 — P.V. y 1 — iirS(y) we obtain the quadratic part of 
the effective action, which contains both real and imaginary parts 
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si%[A,R] = S A [A,R} + S% 
= [ dud 3 k \(AR* + A*R)[- 



(2) 



V 



G(u),k)] + i(A*R - AR*) D(ui, k) + iR*RK(uj,k) 



(14) 



G(uj - 2kv-mv 2 ,k - 2mv) = P.V. / d 3 p 



- mt, 2 , k-2mv)=n J d 3 p S(uj + e f +e % _ f ) [/(/3e ? ) /(/%_ p -) + f(-0e f )f(-/3e^)] 

D{lo - 2kv - mv 2 ,k - 2mv) = ir f d 3 p S(lu + e + H _ f ) [f(Pe f )f(/3e^_ s ) - /(-/3 £?? )/(-/3 e£ _ ff )] 



k—p' 



A term like R*R K(ui, k) with K{uj, k) > in the imaginary part of Sj? is usually resposible for decoherence. It 
suppresses the norm squared of the difference R. In our case we have just a weak inequality K (u>, k) > 0. From its 
definition K > for given u, k iff there exists such a p that u> + ep + £jr_p. = 0. This happens for cj < 2/i — k 2 /4m. 
In the long wave-length limit, for k 2 /m and cj <C /i, this condition is satisfied. The long wave-length modes are 
subject to decoherence and can be considered as classical. 

From now on we restrict our attention to the leading order terms in gradient expansion for v = 



c(2) 
°eff 



d A X {[- 



(2tt 



r dG. 



] + G(0, 0)] (i?*A + c.c. ) + [—(0, 0)] (-iR*d t A + c.c. ) 



V r 1 ^ G mniir/?* d " A 
^ [ 2c^9F (0 ' 0)](jR ; 



c.c.) + [^-(0,0)}(R*d t A 



c.c)+i[K(0,d)]R*R} (15) 



fi.u— 1..3 

As we can see by inverse Fourier transform the propagator of the density matrix contains a factor 

-K(0,S)J dtd 3 xR*R 



(16) 



which is responsible for decoherence in the long wave-length limit. The exponent is minus norm squared of the 
difference R = A + — A_. This term suppresses the off-diagonal elements of the density matrix. The factor ( |l6| ) can 
be rewritten as 



- f dtd^x \ { *1 ~if*R-i£R* 
J L/f(o,o) s s 



The noise has the correlation 



( ?{t, x) e(t', x') ) = K(0, 0) 6(t - t') 5(x - x') 



(2) 

Variation of S e jj with respect to R* leads to 



G 



A 



1 d 2 G 
6 dP 



V 2 A + £ 



(17) 



(18) 



(19) 



where we suppressed the argument (0, 0). It can be checked that < K(0, 0) = 2f3~ 1 {— d u D(0, 0] = ttN/A so that the 
fluctuation-dissipation relation is satisfied; N is density of states at the Fermi level. 

The derivation which leads to Eq. ( |l9| ) is valid above critical temperature, where the expectation value of the order 
parameter is zero. It is based on expansion in powers of A around its vanishing equilibrium value. Below T c the 
coefficient G — (2w) 3 /V is negative and the system is unstable against spinodal decomposition. 



III. BELOW T c 



In the previous section we made an expansion around vanishing order parameter A = 0. Below T c the order 
parameter has a nonzero expectation value which results in a nonzero energy gap for the quasiparticlcs. The gap 
effects can be difficult to generate in perturbative expansion around A = 0. That is why we make a perturbative 
expansion around uniform and time-independent order parameter, A = Ao + 4>, and work out the effective action up 
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to second order in <fi. A first order term vanishes, if Ao satisfies the standard BCS gap equation. The decoherence 
and damping kernels turn out to be nonanalytic at (u),k) = (0,0) and it is not possible to make their gradient 
expansion. Just the nondissipative (reversible) part of the effective equation can be expanded and truncated as a 
partial differential equation; irreversible terms are nonlocal. The limit of T = is considered in more detail in the 
next section. 

The free electronic Lagrangian for a uniform Aq background is 



L„M = / d 



]T (tytM* + ^V&VVa + + (2^) 3 A*^ T + (27r) 3 Ao^]^ 



(20) 



tp'a can be Fourier transformed as in Eq.(0) but a's do not diagonalize the Hamiltonian. We make the Bogolubov 
transformation 

a T (t, k) = u^(t, k) + v^\(t, -k) , 
ai(t, k) = u^i(t, k) - v^(t, -k) , 



w\ 2 - 


i 


1 + _fc 


2 


- e fcJ 


KH 2 = 


f 




2 


- e k\ 




(27t) 3 A 


2u k v k -- 




er 



(2tt)3ASA + e| 



(21) 



is the energy of Bogoliubov quasiparticles of momentum k. It has a gap (27r) 3 / 2 |A |. After the transformation the 
full Lagrangian becomes 

L[ 7 , 4>] = L [j] + L M + Until, 0] , 



LoH - / d 3 k 



lm -- 



d s k 



{illPtla - exilic) 

-(27T) 3 



V 



(A o 0* + Afc 



d 6 k 



^(t,fc)^[ 7 ] + 0(t,fc)F|[7] 



fcL'J ■ ' - ->- k i 

We neglected the the ground state energy terms which depend only on A . The interaction vertex is 
To first order the influence action Sf is 



(22) 



Sp 1 ' = I dtd 3 k 



#.(t,A0(Fg[ 7+ ]) - 0i(*,fc)<fs[7-]> + h.c. 



(23) 



(24) 



where 0^ lives on (anti-)time ordered branch. This first order term cancels with the first order term in L^[cj) + ]— L^[tf>-], 
see Eq.(p2[), provided that A satisfies the BCS gap equation 



2(2tt) 
V 



, tanh 
d A p — 



(25) 



There is no linear term in the effective action if the uniform A is chosen at thermodynamic equilibrium. The 
quadratic term of the effective action is obtained by similar steps as above T c . With the definition <f>± = (f) ± R/2 it 
becomes 
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c(2) _ q(2) q (2) 



dcod 3 k { {<f>R* + (j>*R) 



(2tt) 

V 



G(u,k) 



+ i((f>*R - (f>R*) D(w, k) + iR*R K(u, k) 



[2A^ R*(j>* + /i.e.] U(w, k) + 



(w, k) } . (26) 



The kernels G, D, if, J7, H are listed in Appendix. G and [/ can be expanded in powers of u>, k. D, K, H are nonanalytic 
at (u, k) — (0, 0), they give rise to nonlocal kernels. Any kernel X can be Fourier transformed back 



X[t - t',x~x'} 



dojd 3 k 



e ia,(t-i')+i%-^) X [u,k] 



(27) 



The i?-quadratic terms in e =// can be traded for a complex noise £ and real noises 77^, 77, 

-J d 4 xd 4 x' [R* (x)K(x-x')R(x') + ^(A*) 2 R* (x)H(x-x')R* (x') + ^(A ) 2 R(x)H(x-x')R(x')] 

J DCD^Di h Drj e~I dWx ' [^ K ~^ +r >* (-t^K+'K^H e l / AS+r,i?A ] 



(28) 



Grandient expansion of G, U and subsequent variation of the effective action with respect to R* gives an effective 
equation 



dV iD(x - x 1 ) <j>(x') 



dG 



id t 4>- 



ld 2 G 



6 dk 2 



l d 2 G 
YdJ 2 



0? 



d 2 u 

duj 2 



ld 2 U 



3 dk 2 



2v72 jl* 



A^V 



G 



(27T) 3 



[2f/] Aq + £ + zA* ^ 



(29) 



where we kept derivatives up to second order. The kernels and their derivatives in square brackets are evaluated at 
(oj,k) = (0,0). In the special case of T = 0, which is considered in the next section, one can reconstruct from this 
linearized equation a nonperturbative equation for A(t,x). 



IV. SMALL TEMPERATURE ORT« A | 



The effective equation takes a particularly simple form in the low frequency and long wavelength limit close to 
zero temperature. The nonanalytic kernels (^4|,^5 46) can be ignored provided that temperature is small as compared 
to the quasiparticle gap and \u>\ is less than twice the gap. In this regime there is no decoherence or noise from 
pair-breaking fluctuations. Nontrivial contributions come from the analytic kernels ( |39] , [40| ) . Both of them can be 
gradient expanded. When we restrict ourselves to leading time derivative and potential term, the effective equation 
for small fluctuations will take the form 



id t 4> = 



(2^) 3 4> 



d 3 p 



+ A ^ 



a 



d 3 p 



A o 0*. 







d 3 p 



V ' "7 2e p -[A*A ] ' u ^dA J 2 e # [A£A ] ' u " dA* J 2 e p -[A*A ] ' 
We remind explicitly that e's depend on A . This equation follows from the effective Lagrangian 



L 



eff 



d 3 x F[p] 



id t A 



+ c.c. 



Vefflp] 



def 



The real coefficients arc functions of p — A* A 



= f d 3 p e fi 1_ 

m J 4(2^)3 ep {t, 



x) 



VefM 



(2ir) 3 p(t, x) 
V 



d 3 p 
J2nf 



(30) 



(31) 



(32) 



The coefficients depend on p through the energies e which appear in the definitions ( |39| , ^0| ) but with the energies 
generalized to 
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ep{t, x) = yJ(2n)3p(t,Z) + e} . (33) 
The equation ( |30| ) follows from the Lagrangian (^Tj) as 

* ■ * ' 5L ^ (34) 



where the second variation is around the constant Ao background. The leading phase gradient terms can be recon- 
structed using Galilean symmetry 

L eff = J d?x |-2F[p] (d tX + ^Pj-Vefflp}} > ( 35 ) 

where A == p 1 ! 2 exp i\- This however still leaves undetermined the Galilean invariant term ~ J d 3 x (Vp) 2 . According 
to [jU^D the perturbative version of such a term is — F(po)[\7(5p 1 ^ 2 )] 2 /12mpo, where po = |Ao| 2 . With this term 
included the effective Lagrangian reads 



(36) 



plus higher order derivative terms. 

The Euler-Lagrange which follows from ( |36"1) for a constant A , 6V e ff/5A* = 0, is identical to the zero temperature 
version of the gap equation ( p5| ) . 

We observe that 4F[p] = dV e f / [p]/dp and as such is identified as density of particles. This observation is consistent 
with Ref. ||. In contrast to ^ our Lagrangian ( |35| ) is valid even for nonuniform |A|. 

The effective Lagrangian contains a first order time derivative term ~ dtX- It is n °t a total time derivative. 
This term was neglected in some early treatments because its coefficient F[p] vanishes with the usual approximation 
J d 3 p « 47riV de p , where N is density of states at the Fermi level and wd is the Debye frequency. In ||,|| the 
term is reconstructed with Galilean invariance but only as a total time derivative. This term implies that at zero 
temperature in clean superconductor a moving vortex is feeling the Magnus force transversal to its velocity. 

The gradient terms can not be rearranged as |VA| 2 , there is different stiffness for phase and modulus fluctuations. 
This is in contrast to the situation above T c , compare Eq.(p"9|), and to the Gross-Pitaievskii equation. 

The effective Lagrangian is consistent with the previous work on this topic up to discrepancies which have been 
identified and explained. The new ficture is the p dependence in F(p) and the nonpolynomial character of V e ff(p). 



CONCLUDING REMARKS 



We made a thourough survey of decoherence and noise in neutral s-wave superconductor. The dynamics above 
critical temperature is of Langevin type. In the long wavelength limit the dissipation kernel can be approximated 
by a simple friction term and the source of fluctuations by white gaussian noise. Long wavelength modes decohere 
and are effectively classical. This description survives below T c as long as the gap remains small as compared to 
temperature. As the gap developes with decreasing temperature the dissipation kernel and noise correlations become 
essentially nonlocal - they can not be expanded in powers of derivatives. The decoherence also has a nonlocal kernel. 
As the temperature becomes small as compared to the gap noise, dissipation and decoherence become exponentially 
supressed except for frequencies of at least twice the gap. For low temperatures we obtain a quantum system. 

The quasiparticle gap is the key factor to suppress decoherence and noise. This fact can be easily overlooked if the 
effective theory is constructed by expansion around A = in all temperature regimes. It is difficult to generate gap 
effects in the leading orders of such expansion. 

At near zero temperature we get an effective quantum Lagrangian in terms of the order parameter A. The 
Lagrangian contains a first order time derivative term. This term is responsible for the Magnus force acting on 
a vortex moving with respect to the superfluid and phenomena such as parallel motion of vortex-antivortex pair or 
vortex pair rotation around its center of mass. The existence of this term is also an essential assumption for the 
phenomenological theory of Josephson junctions. 
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APPENDIX 



The nonlocal kernels encountered below T c can be divided into those which have a well defined gradient expansion 
and those which are nonanalytic at (w,k) = (0,0). The analytic kernels are G and U. 



G(cj,k)= P.V. d 3 p{ f([3ep)f(pe~ k _ p ) - f(-{3ep)f(-[3e~ k _ p ) 



/(/%)/(-/%_.) - /(-/3e p -)/(/3e fc -_.) 



\ u v\ 2 \ u k-p\ 2 



w - g p + e fe- P - w + e p ~ e k-p 



As a selfconsistency check we note that for A = the G kernel coincides with the G in Eq. ( ]i~4| ) . 



U(u, k) = ^ P.V. I d 3 p { f{pe p )f{l3e % _ f ) - /(-/?e p -)/(-/?e £ _ p -) 



At zero temperature these kernels simplify to 



e fe-p ' e p 



+ 



e k-p e P 



e P e k-p[( e P + e k-p) 2 ~ «> 2 ] 



G(uj, k) /3 =°° P.V. I d? P 



U(u;,k) =°° P.V.M 



Wp\ 2 \H-p\ 2 



w + ep + eg_ p - u-ep- eg_ # 



d 3 P 



e k-p + e p 



'pH-p^ 2 - ( e P + e k~p) 



The nonanalytic kernels are as follows 

K{u>, k) = TT J d 3 p { [/(/%)/(/%_ # ) + f(-^ep)f(-Pe s _ ff ) 

\ u p~\ 2 \ u k-p\ 2 S ^ + e p + e k-p) + \ v p\ 2 \ v k~p\ 2 <*( w - tp-H-p) 
~f{f}ep)f(-[3e % _ p )+f{-f}ep)f{l3e % _ p ) x 

u p\ 2 \ u k-p\ 2 S ( UJ - e p + e k-p) + \ u p\ 2 \ v k-p\ 2 6 ( UJ + e p~ e k-p) 1 



(37) 



(38) 

(39) 
(40) 



(41) 
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X 



D(u, k) = n I dM[/(/%)/(/%_p) - /(- /3e .)/(- / 3 e ^) 

'f((3ep)f(-(3ej:_ p ) - f(-/3ep)f((3e % _ p ) x 



Again we note that for Ao = these K and D kernels match the corresponding ones in Eq.(14). 



H(u, k) 



7r(27r) 



d 3 P 



e P e k-p 



8(uj + eg + eg_ # ) + S(uj-ep- e £ _ p -) + 



f(pee)f(-/3e s _ p ) + f(-/3e p )f(pez_ f ) 
S(uj- ep + e % _ p ) + S(uj + ep~ eg_ ? )J } 
At zero temperature the nonanalytic kernels become 
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